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CREEP AND CREEP RUPTURE OF STRONGLY REINFORCED 

METALLIC COMPOSITES 


By D.N.Robinson 1 , W.K.Binienda 2 and M.Miti-Kavnma 3 


ABSTRACT: A creep and creep damage theory is presented for metallic composites 
with strong fibers. Application is to reinforced structures in which the fiber 
orientation may vary throughout but a distinct fiber direction can be identified 
locally (local transverse isotropy). The creep deformation model follows earlier work 
and is based on a flow potential function that depends on invariants reflecting stress 
and the material symmetry. As the focus is on the interaction of creep and damage, 
primary creep is ignored. The creep rupture model is an extension of continuum 
damage mechanics and includes an isochronous damage function that depends on 
invariants specifying the local maximum transverse tension and the maximum 
longitudinal shear stress. It is posited that at high temperature and low stress, 
appropriate to engineering practice, these stress components damage the fiber/matrix 
interface through diffusion controlled void growth, eventually causing creep rupture. 
Experiments are outlined for characterizing a composite through creep rupture tests 
under transverse tension and longitudinal shear. Application is made to a thin-walled 
pressure vessel with reinforcing fibers at an arbitrary helical angle. The results 
illustrate the usefulness of the model as a means of achieving optimal designs of 
composite structures where creep and creep rupture are life limiting. 


INTRODUCTION 

When metallic structures operate at temperatures in excess of 0.3 T m , where 
T m is the melting temperature of the metal, it is generally observed that the 
structure suffers time-dependent deformation and eventually ruptures. For components 
made of highly creep resistant materials or where distortion is not a controlling 
factor structural lifetime is often limited by creep rupture. 

In design practice, estimates of creep deformation are often based on 
generalizations of the classical Norton creep law. Predictions of creep failure resulting 
from material deterioration are frequently based on generalizations of the Kachanov 
(1958,1986)/ Rabotnov (1969) damage law. 
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Generalizations of these laws (e.g., Leckie and Hayhurst (1974), Lemaitre and 
Chaboche (1985), Leckie (1986)) have been proposed that account for multiaxial stress 
and that employ damage evolutionary laws reflecting various micromechanisms (Ashby 
and Dyson(1984)). These generalizations commonly employ concepts such as an 
effective stress and an isochronous failure function or surface. Important conclusions 
based on these formulations are that creep deformation, rupture time and failure 
strain all are influenced strongly by the multiaxial state of stress . Further, they 
conclude that damage evolves more or less anisotropically depending on the 
mechanisms involved. 

The objective of this paper is to extend these results concerning creep and 
creep damage to metals that are reinforced by strong unidirectional fibers. The focus 
is on reinforced structures where the fiber direction may vary from point to point 
but that a single material (fiber) direction is identifiable locally (viz., local transverse 
isotropy). The strong fibers are assumed to suppress creep along their direction, but 
in the presence of general states of stress both creep and creep damage occurs. 

The representation of creep deformation follows earlier work concerning strongly 
fiber reinforced metals (Robinson and Duffy (1990), and Binienda and Robinson 
(1990)). Tensorial invariants identified in these works reflecting deviatoric stress 
dependence and the appropriate material symmetry are adopted in the present flow 
law. As the emphasis here is on the interaction of creep and damage ( tertiary creep) 
the relatively small contribution of primary creep (cf.,Petrasek and Titran 
(1988),Grobstein (1989a, 1989b)) is ignored. 

A major assumption here is that the stress dependence of damage evolution is, 
(i) on the maximum tensile stress normal to the local fiber direction and, (ii) on the 
maximum longitudinal shear stress, i.e., the maximum shear stress on planes 
containing the fibers and in a direction along the fibers. It is postulated that the 
effects of these stress components interact to damage the composite internally (along 
the fiber / matrix interface), eventually leading to creep rupture. Under this 
idealization, tensile stress along the fibers causes no creep and no damage; tensile 
stress ( however small) transverse to the fibers causes both creep deformation and 
damage. In reality, stress along the fibers causes damage and leads to a finite 
lifetime, however, the difference in damage rates, and thus rupture times, under 
longitudinal and transverse stress generally exceeds an order of magnitude (Grobstein 
(1989a), Petrasek and Titran (1988)). 
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Studies of creep damage of some monolithic alloys suggest that at lower stress 
levels appropriate to engineering practice, creep rupture occurs mainly by 
intergranular void growth (Cocks and Ashby (1982)). In such cases void growth rates 
at a particular grain boundary are controlled by diffusion mechanisms which, in turn, 
are controlled by the maximum tensile stress transverse to the given grain boundary 
(Hull and Rimmer (1959), Chuang and Rice(1973)). Here, we postulate that the 
fiber-matrix interface in a metallic composite plays a role, on the meso structural 
scale, analogous to that of grain boundaries on the microstructural scale. Thus, at 
high temperature and lower stress levels, we anticipate that the maximum tensile 
stress transverse to the fiber direction may strongly influence void growth at the 
fiber-matrix interface and, consequently may correlate with a creep rupture 
mechanism based on interfacial degradation through diffusion related void growth. It 
is expected that the diffusion mechanisms involved may include those known to 
control void growth rates at grain boundaries as well as other raesomechanisms 
observed to be present in metallic composite systems, e.g., those leading to 
phenomena such as Kirkendall porosity or the evolution of an interphase (Grobstein 
(1989b)). As indicated earlier, account is also taken of the contribution of longitudinal 
shear stress to interfacial degradation. These stress dependencies of damage evolution 
are included by introducing invariants that correspond to the local maximum 
transverse tension and longitudinal shear. These invariants are taken as arguments of 
an isochronous damage Junction (cf.Leckie (1986)). 

As the damage invariants depend on the local stress and material orientation, 
the rate of damage evolution is anisotropic. Thus, rupture time and failure strain 
depend similarly on fiber orientation. The degree of damage is measured by a scalar 
(Kachanov (1958)) and thus does not reflect material directionality. Potentially, this 
may lead to inaccuracies in life prediction of structures subjected to highly 
non-proportional loading; it is expected that this will not seriously affect life 
prediction of proportionally loaded structures, e.g., pressure vessels, (cf.,Ho (1987)). 

Following the development of the creep /creep rupture model, experiments are 
outlined for the determination of the material parameters and functions. Application 
is then made to a simple structure i.e., a thin - walled tubular pressure vessel. 
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THE CREEP/CREEP DAMAGE MODEL 


The form of the creep deformation / damage model proposed by Leckie and 
Hayhurst (1974) is extended for metallic composites. We consider high-temperature, 
isothermal conditions and a composite of given fiber volume percentage. 
Generalization to nonisothermal conditions follows as in Cocks and Ponter (1987). 
Extension to other fiber densities is possible by considering key material parameters 
as functions of fiber volume percentage. 

Invariants identified in Robinson and Duffy (1990) and Binienda and Robinson 
(1990) are adopted to account for material directionality ( local transverse isotropy). 
As earlier, an orientation tensor 

Dij = didj (1) 

is introduced where di (i— 1,2,3) are the components of a unit vector denoting the 
local fiber direction. The invariants Ii and I 2 are defined as (cf., Robinson and 
Duffy (1990)) 

11 = j 2 + -j- 1 2 - j 

1 2 = j- 1 2 

in which 


( 2 ) 


J 2 — — ^ SijSji I — DjjSji 7 — DjjSjkSki (3) 

and sij = the components of the deviatoric stress. The physical meanings of Ii 
and 1 2 are: 


VI7 = S = the maximum transverse shear stress 

VT 2 = <7 = the maximum longitudinal shear stress 
We define an additional invariant as 


where 


and 


T' — < ~ A) + S > 


— <7ii A ~ DijCTji 


r x ; x > 0 

< x > = \ 

1 0 ; x < 0 


( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 
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Physically, 


Jf = the maximum transverse tensile stress (8) 

These invariants enter the model through two functions <j>(crij,Dij) — the flow 
potential function and A(<7ij,Dij) — the isochronous failure function. Guided by 
Binienda and Robinson (1990) we take (j) as: 

* = -§; A-4* w 

where a 0 is a reference stress to be defined subsequently. Consistent with the 
discussion in the previous section, the function A is taken to depend on the 

invariants Jif and df . The surface A( = 1 in the space will be 

termed the isochronous failure surface. A specific functional form of A for a given 
material must be determined by experiment. Possible experiments for this purpose 
might employ reinforced thin-walled tubes of the composite of interest subjected to 
combined tension, torsion and pressure. As pertinent experimental results are not 
available currently, here, in the interest of simplicity, we adopt a linear form (fig.l) 

A = -L ( jf+ a q/) (10) 

in which a is a material constant. 

A modified Calladine and Drucker (1962) potential form serves as the flow law 

i.e., 


Eij - in 9 * 1 

(ii) 

With (j) ,eq. (9), into eq.(ll) there results 


e i j P ij i 

= 2 *»-> 1 
£o °° <P n 

where 

(12) 

Pi] = s ij 2“{3Dij ” ^ij) 

The evolutionary law is taken as (cf., Leckie (1986)): 

(13) 

A" 1 

nAt 0 ' 

in which 

(14) 
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m 


nA , 
X^T ~ 1 


(15) 


The scalar # is the material continuity (Kachanov (1958)); it equals unity for a 
material element entirely intact and zero for an element having lost load carrying 
capacity entirely. 

As indicated earlier, a Q is a reference stress and 
e 0 ,to,A,n,i/ and a 

are material constants; their physical meanings and their determination are the 
subjects of the following section. 

An alternative and useful form of the flow law, eq.(12) can be stated by 
defining an effective stress 


o — -f 2' rijTji 

and an effective inelastic strain 


€ - V ~J~ e ij f ji 

Using eq.(13), the flow law can then be written 


(16) 

(17) 


4 = ( i%r ) = ( 4~ )” ( 18 ) 

Co 

relating the effective stress and effective inelastic strain rate. 

For the limiting case of isotropy , the creep / creep rupture model expressed in 
eqs. (9)-(15) reduces to that proposed by Leckie and Hayhurst (1974) and Leckie 
(1986) for isotropic alloys. This reduction is discussed in the Appendix. 


DETERMINATION OF THE MATERIAL CONSTANTS 

The model expressed in eqs. (9)— (15) predicts no creep and no damage under a 
constant longitudinal tension (fig. 2a). Under a transverse tensile stress 027 — o 
(fig. 2b) , the model reduces to 

4 = (^ r) n < 19 ) 

Co 

* = - ; m = ^r4r - 1 ( 2 °) 
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where e — e 2 2 ~ the transverse creep rate (fig. 2b). With a — cr 0 i,e,. at the 
reference transverse stress, and with = 1 ,i.e., for an undamaged material, eq. 


(19) gives e = e 0 - Thus, the constant e 0 represents the (early) creep rate under 
the transverse reference stress a 0 (fig. 3). 

Again, with a — a 0 eq. (20) can be integrated from $ = 1 to ^ giving 


t A " 1 

* n = (i - 4-)~r 
“0 


( 21 ) 


At failure ( t = tf and ^ = 0 ) we have tf = t 0 from eq. (21). Thus, t 0 

denotes the time to failure under the reference ( transverse ) stress a 0 (fig. 2b). 
Inserting eq. (21) into eq. (19) with a = <r 0 results in 

£ t X " A 

4- = (i - -r)^ - ( 22 ) 


Integrating eq. (22) over t = 0 -* t 0 and e — 0 -» ef , we get for the transverse 
failure strain £f 





— Ato 

(23) 


e 0 



Cf 

(24) 

or 

A = — 


f oto 



The constant A is referred to as the creep damage tolerance in the literature on 
creep of monolithic alloys (cf., Leckie (1986)). In effect, A measures the ability of 
the metal to withstand local straining at points of high strain concentration in a 
structure. In design practice it is desirable that A > 10. 

Now we consider an arbitrary transverse tensile stress a not equal to the 
reference stress a 0 . For ^ = 1 , eq. (19) then reduces to a simple Norton power 


law. With the constants a 0 and e 0 known and data pairs (&,'t) measured from 
the early stages of a set of transverse creep rupture tests, an optimal value of the 
exponent n can be determined. 

Next, we focus on eq. (20) with a not equal to a 0 . Integrating over the 
lifetime ¥ = 1 -» 0 and t = 0 tf we obtain 


li = i 

to (cri'a o) 


v 


(25) 
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Again, with a 0 and t 0 known and the data pairs (o-,tf) measured from a set of 
transverse creep rupture tests, an optimal value of v can be obtained. At this 
point, the only remaining constant to be specified is a in eq. (10). 

Under a constant longitudinal shear stress on — r (fig. 2c), eq. (14) reduces 
to 


Integrating eq. (26) over $ : 


iX-g-)" 1 


a o 

1 -* 0 and t 


_ _ An 
m = j-tt “ 1 

0 -* tf there results 


tf 

to 


{ar/a 0 ) v 


(26) 


(27) 


Combining eq. (25) and eq. (27) and solving for a we get 

t tt 




t ls 

f 


(28) 


where, for purposes of identification, t* fc corresponds to the time to failure under 

the transverse tension a (from eq. (25)) and t* s corresponds to the time to 
failure under the longitudinal shear r (from eq. (27)). If we have several data pairs 

(cr,t^) and (r,t* s ) and v is known, we can find an optimal value of a . In 

particular, if a and r lie on the isochronous failure surface A = 1 then t“ = 

t* s = t 0 and, from eq. (28), a = o/t. This completes the procedure for specifying 
the material constants. 


CREEP RUPTURE UNDER MULTIAXIAL STRESS 


We now examine creep rupture under a constant multiaxial stress. For fixed 
fiber orientation the functions <H <^ij / cr 0 ,D ij ) and A(crij/ cr 0 , Dij) are then constant. 
Eq. (14) can thus be integrated over t = 0 — » t and = 1 -♦ ^ providing 
(cf,eq. (21)) 


= (1 - A 



A— 1 
T" 


(29) 
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( 30 ) 


At failure ( t = tf and $ = 0 ) eq. (29) gives 

tf _ 1 
t 0 A 1 ' 

We recall that for A = 1 , viz., a stress state on the isochronous failure surface, tf 
= t 0 . 

Now consider the alternative form of the flow law expressed in eq. (18). With 
<|> = constant and making use of eq. (29), eq.(18) can be integrated over t = 0 -+ tf 

and 7 = 0 -+ 7f as follows 


tf 

\A 


tf 


1— A 


^ = <H (1 ~ A^4"P~dt 

o to K 0 


(31) 


Using eq. (30) in eq. (31) results in 

7f 


<j> n Atf 


to 


Then using eq. (24), we have 


tf 


tf 


= A1 


(32) 


(33) 


Eqs. (30) and (33) show that the time to failure tf ( crij / cr 0 ,D ij ) and the effective 

failure strain 7f(<7ij/ cr 0 , Djj) are functions of the multiaxial stress and the fiber 

orientation. In principle, for a given stress state the fiber direction can be chosen to 

optimize tf and 7f in accordance with eqs. (30) and (33) and an appropriate 
design criterion. 

In the next section we apply the theory presented here to a fiber reinforced 

thin-walled tubular pressure vessel. There we calculate tf using eq. (30) and 7f 
using eq. (33) for the tubular vessel having a single family of strong helical fibers 
oriented arbitrarily. 
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APPLICATION TO A THIN-WALLED PRESSURE VESSEL 


Consider a thin-walled tube with closed ends that is subjected to an interior 
pressure p. The tube with radius R and wall thickness h is reinforced with a 
family of strong helical fibers making an angle 6 with the circumferential direction 
X 2 (fig. 4). Taking 


( 34 ) 

the stress components a \ j and the deviatoric stress components Sij in the tube 

wall are, respectively 


< 7 ll — ~y~ a 22 = cr 

&ZZ ~ O 12 = O 13 — <723 — 0 


S22 — — S33 — 

S 11 = S 12 = S 13 = S23 == 0 


The components of the orientation tensor Dij are (fig. 4) 

Dn = sin 2 #, D 2 2 = cos 2 #, D 12 = sin#cos# 
D 33 = D 13 = D 23 = 0 

The relevant invariants can then be calculated, 

T _ o 2 
j 2 — 

p 

I 2 — -|-cos 4 # 

giving, from eq. (9) 

♦ = -r 0 J 1 - -h ”* 9 

Additional invariants are 

S = -|-( 2 - cos 2 # ) 

— o ( cos 2 # H — ^-sin 2 # ) 

3 

* = ~r a 

leading to 

— (7 ( 1 — — COS ^0 ) 


(35) 


(36) 


(37) 


(38) 


(39) 


( 40 ) 
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Finally, 


if — — — sin20 


(41) 


and from eq. (10) 

A — (1 g-cos 2 # H — ^-sin20 ) (42) 

Using eqs. (30), (34) and (42) we write the dimensiojiless rupture time vs. fiber angle 
8 as 



(43) 


Evaluating eq. (43) for 8 — 0° (circumferential reinforcement) and 8 = 90° 
(longitudinal reinforcement) and forming the ratio of these results, we get for any 
pressure p, 


t(o°) 

-I = 2 

t(90°) 


(44) 


Thus, the failure time for circumferential fiber strengthening may exceed substantially 

that for longitudinal strengthening. Typically, v s 9 so that 2 U z 512. 

A graph of eq. (43) is shown in fig. 5. The three curves shown are for a — 0, 
1/2 and 1. The first, a — 0, implies that no damage results from the maximum 
longitudinal shear if The last condition, a — 1, implies that if and the 
maximum transverse tension Jf are equally damaging. 

Fig. 5 shows that for a = 0 the rupture time decreases monotonically with 
increasing fiber angle 8. When longitudinal shear if contributes to creep damage, 
viz., a — 1/2 or 1, fig. 5 shows that the rupture time no longer decreases 
monotonically with 8 but exhibits a minimum value at an intermediate fiber angle. 

Combining eqs. (33), (38) and (42) we obtain the dimensionless effective failure 
strain vs. 8 , viz., 



(J 1 --j-cos^) n 
— 2 — cos 2 $ + — ^ — s i n2 0 } 


(45) 
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Similarly evaluating eq. (45) for 9 = 0° and 6 
results, we have for any pressure p. 

7f(0°) 

_ 2^~ n 

ff(90°) 


90° and taking the ratio of these 


(46) 


Now for definiteness we consider specific values of the exponents n = 6 and v — 
9 and take the pressure as pR/cr 0 h = 1. The result is plotted in fig. 8. As earlier, 
three curves are shown corresponding to a = 0, 1/2 and 1. With the given values 

of p,n and v each curve shows 7f = e f at 6 = 90° i.e., the effective failure 

strain is equal to that corresponding to transverse tension at the reference stress a 0 . 

At 9 = 0° each curve shows ef = 2 9 " 6 ef = 8 €f, consistent with eq. (46). 

For the curve a — 0 in fig. 6, 7f reaches a maximum of al4 tf at 

9 s 25°; conversely, the curves a — 1/2 and 1 each exhibit a minimum at an 

intermediate fiber angle. 

The simple structure considered in this section is not of essential practical 

importance but serves as an illustrative example of the dependence of tf and 7f 

(or equivalently, an effective creep damage tolerance X = (7f/tf) A ) on the 
multiaxial ( biaxial) stress state and on the fiber orientation. In a general design 

situation we may want to maximize both tf and 7f (or X), i.e., we may want a 
long structural lifetime and a large tolerance to local straining at possible strain 
concentrations. In our present example, we see by comparing figs. 5 and 6 that these 
maxima do not always occur at a common fiber angle 9 , suggesting the necessity 
of a design compromise in the selection of the optimal fiber orientation. It is 
intended that the simple model presented here may provide guidance for the 
practicing structural engineer in making such design considerations. 


SUMMARY AND CONCLUSIONS 

A creep and creep damage model is presented for metallic composites with 
strong fibers. The intended applications are reinforced structures in which the fiber 
direction may vary throughout but a single fiber direction can be identified locally 
{local transverse isotropy). 
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The representation of creep deformation follows earlier work and includes a flow 
potential function <j) that depends on invariants reflecting the local stress and 
material orientation. As the focus is on the interaction of creep and damage, the 
relatively minor role of primary creep is ignored. 

The creep rupture model is a generalization of classical continuum damage 
mechanics and includes an isochronous damage junction A taken to depend on 
invariants specifying the maximum tensile stress normal to the local fiber direction 
and the maximum longitudinal shear stress. These stress components are assumed to 
damage the composite at the fiber/matrix interface, eventually causing creep rupture. 
It is conjectured that at high temperature and low stress the fiber/matrix interface 
plays a role, on the meso structural scale, analogous to grain boundaries on the 
microstructural scale. An anisotropic damage mechanism is posited where interfacial 
degradation occurs through diffusion ( and hence stress) controlled void growth, e.g., 
Kirkendall void growth. 

A procedure is outlined for the determination of the several material parameters 
on the basis of creep rupture testing under transverse tension. Additional testing is 
required to determine the functional dependence and form of the isochronous damage 
Junction. Here, a linear form is assumed requiring a limited number of creep rupture 
tests under longitudinal shear. 

General results are obtained for the time to failure tf and the final effective 

failure strain If under constant multiaxial stress; each is found to depend not only 
on stress, but on fiber orientation. Thus, for a given stress state, the model permits 

the selection of, a fiber orientation that may optimize tf and/or tf according to a 
relevant design criterion. 

The theory is applied to a thin-walled tubular pressure vessel with strong fibers 

at an arbitrary helical angle. The quantities tf and If are calculated for the thin 
tube as functions of fiber angle. As anticipated, the rupture time tf is a 
maximimum for circumferential reinforcement of the tube. Depending on the relative 
strength of the contributions of transverse tension and longitudinal shear to the 
evolution of damage (viz., the form of the isochronous damage function A ) tf 
may reach a minimum at an intermediate fiber angle (fig. 5). The effective final 

strain ef can achieve a maximum or a minimum at an intermediate angle, again 
depending on the relative importance of transverse tensile and longitudinal shear 
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damage (fig. 6). In general, simultaneous optimization of tf and If results in a 
design compromise. 

The example problem illustrates the potential usefulness of the proposed model 
as a means of realizing optimal designs of composite structures where creep and creep 
rupture are life limiting. Of course, the accuracy of predictions made using the theory 
depends intimately on the forms of the functions <j> and A and on the material 
constants. Correct choices of these functions and parameters for a given composite 
material rest on careful experimentation. Some experimental justification for the 
choice of <|) is provided in Binienda and Robinson (1990). However, the functional 
dependence and form of the isochronous damage function A is largely conjectural, 
particularly the linear form adopted here, and must be verified on the basis of 
detailed phenomenological and microstructural observation. Key verification 
experiments are currently being defined. 
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APPENDIX B. ISOTROPIC LIMIT 


The isotropic limit of the model expressed in eqs. (9) - (15) is found by taking 

Djj to be an isotropic tensor with Dh = 1 , i.e. , Dij = -g-£ij. Under this limit 
<j) in eq. (9) becomes 


and Tij in eq. (13) is 


<j) - -J-V J2 


Tij — Sij 


(IB) 


(2B) 


Following Binienda and Robinson (1990), we define a new reference stress a 0 ' and a 
corresponding creep rate e 0 ' as 

(To' = 


e 0 ' 


2 • 
fo 


•vT3" 


(3B) 

(4B) 


With these definitions <() ,eq. (IB), becomes 


d> = iL$h 

and the flow law eq. (12) reduces to 


(SB) 


c ij 


Co' 


3 d>"-‘ Sij 1 
~T <t> 


O o' 


ipn 


(6B) 


as in Leckie and Hayhurst (1974) and Leckie (1986). 

Similarly, eqs. (14) and (15) reduce to those of Leckie (1986) with the function 
A ,eq. (10), taken as 

O’max 


a = 4 = 05 


or 


A = 


<V 


(7B) 


where <7 max = the maximum (tensile) principal stress. 

We note from eq. (3B) that the ratio of the reference stress cr 0 ' of the 
isotropic limit to the ( transverse ) reference stress a 0 of the anisotropic composite is 

v/ 3 / 2. This result is expected in that the material deforms as a J 2 (von Mises ) 
material in the isotropic limit (cf.,eq. (5B)) whereas under the constraint of the 
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strong fibers it is forced to deform in transverse tension as a maximum shear stress 
( Tresca ) material. The creep (or yield) strengths of these material idealizations are in 

the same ratio ^ 3 / 2. 
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Figure 2. Stress states: (a) longitudinal tension, (b) transverse tension, 

(c) longitudinal shear. 
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Figure 



o o 


Figure 3. Creep response under reference transverse stress a 0 . 



Helically reinforced thin-walled tube under interior pressure p. 
local fiber direction. 


d denotes 
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Figure 5. Dimensionless rupture time vs. fiber angle 0 . 
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